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We give the details of the analysis for critical properties of spin-gap phases in one-dimensional lattice 
electron models. In the Tomonaga-Luttinger (TL) liquid theory, the spin-gap instability occurs when 
the backward scattering changes from repulsive to attractive. This transition point is shown to be 
equivalent to that of the level-crossing of the singlet and the triplet excitation spectra, using the 
c = 1 conformal field theory and the renormalization group. Based on this notion, the transition 
point between the TL liquid and the spin-gap phases can be determined with high-accuracy from 
the numerical data of finite-size clusters. We also discuss the boundary conditions and discrete 
symmetries to extract these excitation spectra. This technique is applied to the extended Hubbard 
model, the t-J model, and the t-J-J' model, and their phase diagrams are obtained. We also discuss 
the relation between our results and analytical solutions in weak-coupling and low-density limits. 
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I. INTRODUCTION 

Spin-gap phases of one-dimensional (ID) electron sys- 
tems have been studied for long times. This research has 
been motivated by the phase transitions in ID organic 
conductors. In the past decade, the discovery of high-T c 
superconductivity strongly stimulated this study. 

The spin-gap transition in ID lattice models had been 
mainly analyzed by two approaches: The one is the weak 
coupling theory based on the bosonization theory and 
renormalization group. The other is numerical calcula- 
tion in finite-size systems which is free from approxima- 
tion. In the former scheme, the existence of the gap is 
argued by investigating the backward scattering effect on 
the fixed point, but the validity of the result is ensured 
only in the weak coupling limit. On the other hand, in 
numerical calculation, the analysis is done by a direct 
evaluation of the gap and the finite-size scaling method. 
In this approach, a singular behavior of the gap near the 
critical point makes it difficult to make out the instabil- 
ity. 

In order to illustrate the difficulty in the determina- 
tion of the phase boundary, let us consider the Hubbard 
model, 



Hum = -ty2(c\ s Ci + i iS + H.c) 



(1) 



This model has a spin gap for U < 0. According ta 
the Bethe-ansatz result for the charge gap at half- filliiigEJ 
combining a canonical transformationo, we can obtain 
the asymptotic behavior of the spin gap near the critical 
point U = as 



Since the gap opens slowly near the critical point, it is 
very difficult to find the critical point using conventional 
finite-size scaling method. 

In this paper, we give a remedy for this problen 
The many-body problem is often simplified by using the 
notion of universality. Generally, ID electron systems 
belong to thc^iiyprsality class of Tomonaga-Luttinger 
(TL) liquidsB'Sfl'tffl which are characterized by gapless 
charge and spin excitations and power-law decay of cor- 
relation functions. This behavior can be described by 
the bosonization theory or the c = 1 conformal field the- 
ory (CFT). In this scheme, the phase transition to the 
spin-gap phase is understood as an instability caused by 
the backward scattering process using the renormaliza- 
tion group technique^!, and a spin gap opens when the 
backward scattering turns from repulsive to attractive. 
This transition point is equivalent to-tksrlevel-crossing of 
the singlet-triplet excitation spectratlHl3ll3, by taking ac- 
count of the logarithmic corrections originated from the 
backward scattering. 

In this paper, we will analyze the following models 
based on this notion. The first example is the extended 
Hubbard model which is given by 



^EHM = Hhm + V ^ 



mn i+ \. 



(3) 



AE ~ y/2t\U\e 



For the study of spin-gap transitions, this model has been 
analyzed by the g-ology for weak coupling regionM. The 
numerical calculation was performed by tbergxact diag- 
onalization with finite-size scaling methodoli§t!l How- 
ever, the spin-gap phase boundary has not been clarified. 
The next example is the t-J model described by 



TTt/\U\ 



(2) 



1 



H t -j = -f^(ct s c i+ i, s +H.c.) 

is 

+jJ2( s * ■ - mm+i/A), (4) 



quantum numbers. The second is derivation for the log- 
arithmic corrections. The third explains the calculation 
in two-electron systems. 



where Ci S = c; s (l— n^-s). This model is obtained by dop- 
ing holes in the Heisenberg spin chain. For this model, 
the weak coupling treatment is difficult due to this strong 
coupling constraint, however, the universality class of 
this model is known as TL liquids, from the analysis for 
the exactly solvable cases at J/t = 0|-(p4nless fermion) 
and J/t — 2 (super-symmetric point )EZHla. The remain- 
ing region was analyzed using the exact diagonalization 
by Ogata et alS3. Their phase diagram shows the en- 
hancement of the superconducting correlation (K p > 1) 
and the phase separation (K p — > oo) for the large J/t 
region. According to their result, the spin-gap phase 
does not exist except for the low density region. Vari- 
ational approaches also played roles in the analysis for 
this modelEjl£j'E3, but could not establish clear solution 
for the spin-gap phase. 

Extensions of ±ip_t^Ljiiodfil are also considered by 
many researcherscjlHJ'OLjO'La. In spin systems, a spin 
gap opens by the effect of frustration or dimerization. 
Metallic spin-gap phases can be generated by doping 
holes in these spin systems. In this paper, we concen- 
trate our attention on . the J- J- J' model which includes 
the effect of frustratior£3Eil: 



H 



t-j-j' 



— Ti-t-j + J'^^i^i ■ Si+2 



n t n. l+2 /4:). (5) 



We introduce a parameter a for the strength of the frus- 
tration given by a = J' /J. At half-filling (n = 1), this 
model becomes an S = 1/2 frustrated spin chain. In this 
case, the ground state at a = 1/2 is the two-fold degener- 
ate dimer state with ajsoin.gap, and the ground state en- 
ergy density is — 3/4J^'E2rEil. The fluid-dimer transition 
occurs at a c — 0.2411E3. Upon doping of holes, the sys^. 
tern may become metallic, and the spin gap is reducedO 
but persists for the finite doping. The phase diagram of 
this model for n ^ 1 at a — 1/2, using the exact diago- 
nalization, was obtained by Ogata, Luchini, and RiceEJ, 
but the phase boundary of the spin-gap phase was also 
remained to be ambiguous. 

This paper is organized as follows. In Sec.|J, we 
discuss, based on the continuum field theory, that the 
level crossing of singlet and triplet excitation spectra 
gives the critical point of the spin-gap transition. In 
Sec. Ill, we consider boundary conditions for the unique 



ground state, and discrete symmetries of wave functions 
to identify the energy spectra observed in our analysis. 
In Sec. IV, we analyze representative models introduced 
above, and clarify the spin-gap region in the phase dia- 
grams, and check the consistency of our argument. Fi- 
nally, in Secjy|, we present our conclusions. 

The paper also contains three Appendices. The first 
shows the relation among the different notations for the 



II. CONTINUUM FIELD THEORY 



A. Effective Hamiltonian 



Let us start our argu 



m the Abelian bosoniza- 
The low-energy excita- 



tion theory of electronsBlI^t 
tions are described by continuous fermion fields which 
are defined by 



Cj, 3 -* ll>-L,s{%) + i>R,s( X ) ( 6 ) 

The boson representation of the fermion operator is 
1 



_ e irfeFXgi/\/2-[r(0p+s0 <T ) — 6 P — s8„] 



(7) 



where a is a short-distance cutoff, r — R, L and s =|, | 
in that order. The phase fields are defined 



refer to 




as 




(j) v (x) = 


m 




e v {x) = 


ITT 





^2 A p( x ) Kb) + ^lO)] - — ^- h v 



L 



(8a) 



M x ) Kb) ~ Mp)] + —j—m v , (8b) 



p^O 



L 



where A p (x) = ^e 1q IpI/ 2 'p 31 , and v r is the charge (y = 
p) or the spin (y = a) density operator. These phase 
fields satisfy the relation \<p v {x), 9 v (x')] = — i7rsign(x — 
x')/2. 

Using above relations, effective Hamiltonian of a ID 
electron system is described by the U(l) Gaussian model 
(charge part) and the SU(2) sine-Gordon model (spin 
part), 



TL = Tip + TLo 



(2na) 



dx cos 



0- 



(9) 



Here gi± is the backward scattering amplitude and for 
v = p,<J 



^ = t 1 / dx [ K »(W 2 + K-\d x d> u f} 

Zi: Jo 



(10) 



where v v and K v are the velocity and the Gaussian cou- 
pling, respectively, for the charge (y = p) and the spin 
(y = er) sectors. In the TL phase (gi± > 0), the pa- 
rameters K a and are renormalized as K* — 1 and 
g 1± = 0, reflecting the SU(2) symmetry. 

The phase fields defined in eqs.(||) satisfy the following 
boundary conditions, 



(j> y (x + L) = cj>v(x) — V2irn u , 
6 v (x + L) = v {x) + V2-Km v . 



(11a) 
(lib) 
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The quantum numbers m„ and n v are defined by the 
eigen values of the total number operators iV ri g (mea- 
sured with respect to the ground state) for right and left 
going particles (r = R, L) of spin s 



[(N m + N^)±(N Rl +N u )]/2, (12a) 
[(N R1 -N L1 )±(N RL -N LL )}/2. (12b) 



Here the upper and lower sign refer to charge (y = p) 
and spin (y = a) degrees of freedoms, respectively. Thus 
n v denotes excitations involving the variation of parti- 
cles numbers and m„ indicate current excitations. If we 
require N r>s to be an integer, the possible value of the 
quantum numbers are restricted as 

This is the selection rule for the quantum numbersBEl. 

B. Excitation Spectra and Boundary Conditions 

First, we consider the excitation spectra for gi± = 
case. If the system is periodic, and has unique ground 
state, the ground state energy of the system with length 
L is given byO 



J^o{L) - Le Q — c, 

6L 



(14) 



where the central charge c characterizes the universality 
class of the model. The finite-size corrections for the exci- 
tation energy and momentum of the system are described 
byffl 



2-KVp 2-KV a 

& — J^n = 



P - Pq = j-(s p + s a ) + 2m p k F , 



(15) 
(16) 



where kp = -kN/2L is the Fermi wave number. x v = 
A+ + Aj7 , s u = A+ — A~ are the scaling dimension and 
the conformal spin, respectively, where the conformal 
weights for each sector are given by 




+ < 



(17) 



Here the integer ri^ denote descendant fields which de- 
scribe particle-hole excitations near the Fermi points. 
The scaling dimensions are related to the critical expo- 
nents for the correlation functions as 



(18) 



Therefore, there is one to one correspondence between 
the excitation spectra and the operators. The operators 
correspond to the excited states are given by 



O 



rn .m CT ,n c ,n c 



cx e 



iv / 2(in p ()ip-|-m o -i;l er -|-npep+ra tT 0„) 



or their linear combinations. From eqs. (0) and @, the 
Fermi operator takes the following boundary conditions 
depending on the excited states: 



ip riS (x +L) = ip riS (x)j^ m >> +m <' +n >> +n ° ) . 



(20) 



This means that the excited states given by arbitrary 
combination of quantum numbers are realized by chang- 
ing the boundary conditions, while, for fixed boundary 
conditions, the possible excited states are restricted by 
the selection rule (|l5|). 

The excitation spectra on which we will turn our at- 
tention can be obtained based on the operators for the 
charge-density-wave (CDW) and the spin-density-wave 
(SDW): 

CCDW = ^Ii^RT + V'LlV'Rl 

= — exp(i2fcpa; + iV2<fi P ) cos(v / 2</> (7 ), (21a) 
na 

Osuw.z = i/it^RT - K|^EU 

= — exp(i2fc F a; + iy/2<j> p ) sin(V20 CT ), (21b) 
na 

CSDW,+ = Ipl^Rl 

= — !— exp(i2fc F a; + iV26 p ) exp(+iV20 a ). (21c) 
2na 

These excitation spectra consist of the charge part 
which carries the momentum 2fcp, and the spin 
part which forms singlet (v^cos-y/^o-) and triplet 
(V2 sin V24> a , exp(±iv / 2^ CT )) states. Note that the spin 
part of the singlet and the triplet superconducting oper- 
ators (SS, TS) are obtained with fcp = and replacing 
cj) p -> 6 p . 

If charge-spin separation occurs, the spin excitations 
in eqs.([2l|) (m CT = 1 or n„ = 1, otherwise = 0) can be 
extracted by using anti-periodic boundary conditions fol- 
lowing eq. (|2C|) ■ In the continuum field theory based on 
the TL model, the dispersion relation is approximated 
by linearized one, so that the deviation from the approx- 
imated dispersion become smaller if the excitation ener- 
gies become lower by eliminating the charge excitations. 
Therefore, the precision of the analysis is enhanced by 
twisting the boundary conditions. 

The twisted boundary conditions are also important in 
identification of excitation spectra. Under anti-periodic 
boundary conditions, the momenta of these states are 
reduced to zero. Then we can define the parity trans- 
formation to classify these spectra. Although the space- 
inversion operator and translation operator do not com- 
mute, we can classify these spectra simultaneously by 
wave numbers and parities, if the wave number k takes 
or 7r. From eq. (]?]), the phase fields </>„ change under the 
parity (V: R^-L), and the spin-reversal transformations 
(T: T~!) 



(19) 



V 
T 



(22a) 
(22b) 
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Thus operators have discrete symmetries as V = T = 1 
for the singlet (v2cos and V = T = — 1 for the 
triplet with S z = (y/2 sin \/2(f> a ). The discrete sym- 
metries of the wave functions of these excited states are 
determined by combinations of those of the ground state 
and the operators. Further discussion for the discrete 
symmetries will be given in the next section. 



C. Renormalization Group 

Next, we consider the renormalization (gi± 7^ 0). By 
the change of the cut off a — > e dl ^ the coupling constant 
gi± and K a are renormalized as 



dj/o(0 
dl 

dl 



-yo(i)yi(i), 



(23a) 
(23b) 



where y (0 = 2(1^^ - l),yi(l) = gi±./itv a . For the 
SU(2) symmetric case yo(l) — yi{l) (the level- L-SJJI2) 
Wess-Zumino-Novikov-Witten (WZNW) modeS&EI) 
and ?/o (0 > 0, the scaling dimensions of the operators 
for singlet and triplet excitations splitJogarithmically by 
the marginally irrelevant coupling aaHJO (see Appendix 



B) 



riet 



triplet 



fJO 



1 



2 Ay Q \nL + l 



1 1 



fJo 



2 4y lni + l 



(24a) 
(24b) 



where y = yo(0) and we have set I = InL. When y < 0, 
yo(l) is renormalized to yoQ) —> — oo, then a spin gap 
appears. At the critical point (yo = 0), there are no 
logarithmic corrections in the excitation gaps (the loga- 
rithmic correction from higher order also vanish) . There- 
fore, the critical point is obtained from the intersection 
of the singlet and the triplet excitation spectrallini3'tj. In 
this casep— sve can determine the critical point with high 
precisioriij, sinoe-ithje remaining correction is only x v = 4 
irrelevant fieldscilca. This irrelevant field, which does 
not exist in the pure sine-Gordon model, comes from the 
nonlinear term neglected when linearizing the dispersion 
relation near the Fermi level in the course of the bosoniza- 
tion. 

The physical meaning of this transition point (jjo — 
0) is the one where the backward scattering coupling 
changes from repulsive to attractive. Moreover, at the 
critical point, the SU(2) symmetry is enhanced to the 
chiral SU(2)xSU(2) symmetry^, since the spin degrees 
of freedom of the right and the left Fermi points become 
independent. 

Eq.@ also explains the fact that the SDW (CDW) 
correlation is dominant for K p < 1 region with(out) spin 
gap, while for K p > 1, the TS (SS) correlation is domi- 
capCJ 



nant with (out) spin gap 



Finally, let us consider the massive region. The be- 
havior of the gap is explained from the two-loop renor- 
malizatjanpgroup equation of the level-1 SU(2) WZNW 
modeB&SEl 



dyo(Q 
dl 



1 



VW - o2Af (0 



(25) 



If we define the correlation length £ as yo(hi£) = — 1 and 
the energy gap as AE — v a /^, then one can derive the 
asymptotic form of the spin gap by solving the differential 
equation for |yo(0l ^ 1 as 



AE (x 



-Const. / 1 2/o | ) • 



(26) 



Note that eq. ( p6| ) is the same asymptotic behavior as the 
spin gap of the negative-?/ Hubbard model at half-filling 
given by eq.(§). 



III. UNIQUENESS OF GROUND STATE AND 
DISCRETE SYMMETRIES 

In the previous section, the ground state is assumed to 
be a singlet, so that we should consider the way to make 
the singlet ground state in the finite-size systems. Fur- 
thermore, we also discuss the discrete syiptnetries of wave 
functions to identify the energy spectrao. The symme- 
tries depend on the choice of representations for wave 
functions, so that we consider in the representative two 
cases: one is the standard electron systems such as the 
(extended) Hubbard model. The other is doped spin 
chains like the t-J model. In the following argument, 
the electron hopping is restricted to the nearest neigh- 
bor, and the number of electrons is assumed to be even. 
The results are summarized in Table [|. 



A. Hubbard-type Models 

It is convenient to use the following representation of 
the basis to describe the Hubbard-type models which per- 
mits double occupancy: 

|*A> = /a ("4, ■ • • ,n M ;n M +i, ■ ■ ■ ,n N ) 

n 1 <---<n M --n M + l<---<n N 
M N 

i=l j=M+l 

where 1 < m < L and the periodicity m + L —>■ rii is 
assumed. N is number of electrons, and M is number of 
electrons with down spins. We call this representation as 
"basis A" . In the case of the (extended) Hubbard model, 
all off-diagonal matrix elements of the Hamiltonian arise 
from the hopping term. Then the sign of these elements 
are negative as far as no electron hops across the bound- 
ary. When an electron moves across the boundary, the 
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sign of the hopping amplitude changes depending on M 
reflecting anti-commutation relations of the Fermi oper- 
ators. If the periodic boundary conditions are assumed, 
the hopping amplitude at the boundary become +t for 
M = even case, and —t for M — odd case. ._. 

According to the Lieb-Schultz-Mattis theoremcll (the 
Perron- Frobenius theorem), if all off-diagonal elements 
of a real symmetric irreducible matrix are non-positive, 
then the all vector elements for the lowest eigen value 
have the same sign. Therefore if the lowest eigen state is 
degenerate, these states can not be orthogonal. Thus the 
ground state is proved to be a singlet. In order to real- 
ize this situation in the (extended) Hubbard model with 
basis A, we choose anti-periodic boundary conditions for 
M = even case, and periodic boundary conditions for 
M = odd case. Then all off-diagonal matrix elements 
become non-positive and the ground state is proved to 
be a singlet. This selection of the boundary conditions 
are equivalent to those derived from the Bcthc-ansatz eg.- 
sult of the Hubbard model in the strong-coupling limito. 

The t-J model (0) can also be described by using 
the basis A. This Hamiltonian has off-diagonal matrix 
elements originated from the exchange interaction, in 
addition to the hopping term. The exchange process 
gives — J/2, where the negative sign arises from the anti- 
commutation relations of fermions. Thus the all off- 
diagonal matrix elements are non-positive if the bound- 
ary conditions are chosen as discussed above. This sit- 
uation does not change if the three-site term is added. 
Therefore the ground state of the t-J model is also proved 
to be a singlet. 

Now we define an operation for the parity transforma- 
tion (space inversion) as 

V 7a ("ir • • ,n M ;n M+ x,- ■ ■ ,n N ) = 

/a("m, • • • , fii;n N , ■ ■ ■ ,n M +i), (28) 

where fij = L + 1 — rij. The spin-reversal transformation 
can be defined only for M = N/2 cases as 

Tfafai,- ■ ■ , n M ;n M +i, ■ ■ ■ ,n N ) = 

/a(«m+i, • • • , n N ; ni, ■ ■ ■ , »Af). (29) 

The eigen values of the operator V and T can take ±1. 
Since the all vector elements of the ground-state wave 
function have the same sign, the ground state discussed 
above satisfies V = T = 1 and k = 0. 

The discrete symmetries of wave functions for the ex- 
cited states are determined by combinations of those of 
the ground state and the operators given by the bosoniza- 
tion argument. Thus V = T = 1 for the singlet 
(v^cosv^^ct), and V = T = — 1 for the triplet with 
S z = (\/2 sin \f2~4> a ). On the contrary, for the triplet 
state with S z = ±1 (exp(±iv20er))i the symmetries are 
the same as those of the ground state (V = 1, k = 0), 
because the all off-diagonal matrix elements are non- 
positive. In this case, the SU(2) symmetry seems to be 
broken, this discrepancy of the parity in the triplet states 



is due to the definition of eq. (|28|) which does not imply 
the change of the sign stems from the permutation of 
fermions in the V ,T transformations. If we take account 
of the anti-commutation relation of the Fermi operator 
in these transformations, then we get V' = (— 1) M V and 
T = {-1) N / 2 T. Thus the SU(2) symmetry is recovered. 

B. Doped Spin Chains 

On the other hand, the models with a constraint to 
eliminate doubly occupied sites, such as the t-J(-J') 
model, are obtained by doping holes into S — 1/2 spin 
chains. These models are described by the "basis B" 
which is defined by 

N 

I*b) = X! fB{ni,s 1 ;---;n N ,s N )'[[cl HSi \va,c). 

tli<"'<Iljf i=l 

(30) 

Hereafter, we argue based on the t-J model. In order 
to make the off-diagonal matrix elements stem from the 
exchange process non-positive, we introduce a new basis 
/' with a sign factor as 

/ b (...)=(-1)^^' + 1/2) /b(---): (31) 

where Sj — ±1/2. Note that j denotes the coordinate of 
the squeezed spin systemc3. If an electron moves across 
the boundary, this sign factor changes as 

(_l)Er=i^^+V2) ^ i ('' +1/a >(-l) af . (32) 

Therefore, an additional negative sign appears if M — 
odd. In addition to this, a negative sign is also added, 
reflecting anti-commutation relations of Fermi operators, 
so that the all off-diagonal matrix elements become non- 
positive if we chose anti-periodic boundary conditions for 
M = even, and periodic boundary conditions for M = 
odd. Then all /g have the same sign so that the ground 
state is proved to be a singlet. In this condition, the 
sign of /b does not change by the shift operation (rij — > 
rij + 1), so that the wave number of the ground state is 
k = for any M. 

Next, we define parity and spin-reversal transforma- 
tions as follows: 

V f-a(nx, sx\ • • • ; n N , s N ) = fn(nN, si; • •• ; fii,sjv), (33) 
T/g(m, si; • • • ; tin , sn) = /b(«i, — si; • • • ; hat, -sat), (34) 

where fij = L+l — rij, Sj = s N+i—j ■ As we have proved , 
all / B have the same sign in the ground state, so that we 
only consider the variation of the sign factor of eq.(|3l|) in 
these transformations. One can easily show that the par- 
ity and the spin-reversal transformations bring additional 
negative sign only for M — odd case. Therefore, the sym- 
metries of the ground-state wave function are V = T = 1 
for N/2 = even, V = T = -1 for N/2 = odd. 
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The symmetries of the excited states can also be dis- 
cussed in the same way as in the case of the basis A. In the 
basis B, V = T = ±1 for the singlet (\/2 cos \/2<j) a ) and 

V = T = qpl for the triplet with S z = (y/2 sin \/2><r), 
where the upper (lower) sign denotes N/2 — even (odd) 
case. The triplet states with S z — ±1 are V = qpl 
(exp(±iv / 2^<r)). Note that SU(2) symmetry is conserved 
in the parity of triplet excitations due to the M depen- 
dence of the parity. 

For half-filling (N — L), the ID t-J model should be 
equivalent to the 5=1/2 Heisenberg spin chain. The 
discrete symmetries known for the Heisenberg chain are 

V = T = 1, k = for L/2 = even and V = T = -1, k = 
7r for L/2 = odd. In this case, the negative sign in cq.(^) 
for M — odd case can be canceled by wave number k = n, 
instead of changing the boundary conditions. Therefore 
the symmetries in these two cases are consistent. 

Unfortunately, the above proof can not be applied to 
the t-J- J' model which includes the anti- ferromagnetic 
next-nearest-neighbor interactions. However, boundary 
conditions and discrete symmetries of this model are ex- 
pected to be the same as those of the t-J model, as far 
as no instability takes place. 



and the Drude weight is given by the excitation m p 
<i>/7r as 



D 



K P Vp 



(38) 



Therefore, K p is obtained as K p = n\J Dn 2 n/2. 

In order to obtain scaling dimensions of the spin de- 
grees of freedom, the spin-wave velocity is calculated by 
the following relation, 

E(L,N,S=l,q = 2n/L)-E (L,N) 
V '=£So 2ttJl ' (39) 

The extrapolation is done by the function v a (L) — 
v a (oo) + A/L 2 + B/L 4 . These corrections are explained 
by the irrelevant fields with x v = 4. 

In the following, we analyze some models. Since 
there are too many instabilities in the extended Hub- 
bard model, we consider the t-J model first to turn our 
attention on the spin-gap instability. 



A. The t-J Model 



IV. PHASE DIAGRAMS OF THE LATTICE 
MODELS 



Now we start our analysis for the models introduced 
in Sec|| Besides the spin-gap instability, the charge de- 
grees of freedom is described by the single parameter K p . 
For K p > 1, the superconducting correlations dominant, 
while for K p — > oo, a phase separation takes place. In or- 
der to determine K p , we need two independent physical 
quantities. In our case, we calculate the compressibility 
and the Drude weight--. In finite-size systems, the com- 
pressibility is given byc3 

L ( E (L, N + 2) + E (L, N - 2) - 2E (L, N) 

(35) 

where E (L,N) is the ground state energy of a system 
with size L and N electrons (n = N/L). On the other 
hand, the Drude weight is given by the relatior 



L d 2 E (<S>) 

2 



(36) 



where $ is the flux which penetrates the ring. In the 
continuum field theory, these two physical quajaiitics 
are described by the parameters of TL liquidsoEj'EJ. 
The compressibility is given as excitation n p = ±1 in 
eqs.©,©: 



1 



2K n 



(37) 



Here we analyze the spin-gap phase of the ID t-J 
model (||) by the above explained method. We diago- 
nalize L = 8-30 systems by the use of the Lanczos and 
the Householder algorism. 

In order to investigate the structure of excitation spec- 
tra in detail, we show in Fig.[| the spectral flow (flux de- 
pendence of energy) of N/L = 4/8 system at J/t — 2. In 
this case, the boundary conditions are fixed to the ground 
state, so that the singlet and the triplet excitation spec- 
tra appear at $ = n which is equivalent to the twisted 
boundary conditions. However, the wave number is not 
k = but k = 2fcp. This momentum-shift is explained 
by the relation fc($) = k(Q) + N&/LB. The singlet and 
the triplet excitation spectra are connected adiabatically 
from those of the CDW and the SDW, respectively. If the 
linearized dispersion relation is exact, these two spectra 
move parallelly versus <& in this diagram. As shown in 
eqs.(|2l|), the charge degrees of freedoms contribute the 
same amount in the CDW and the SDW spectra, so that 
the critical point can be obtained by the level crossing 
of these spectra. However, in this case, the precision 
become lower due to the irrelevant fields, and the iden- 
tification by the parity become impossible due to the in- 
commensurate wave number. 

Figure shows the singlet and the triplet excitation 
spectra (L = 16, n = 1/2) versus J/t. The level crossing 
takes place at J/t ~ 2.7. The size dependence of the 
critical point is shown in Fig.|[ Since the critical point 
is almost independent of the system size, the phase di- 
agram can be constructed without extrapolation. Then 
we obtain Fig.^J(a). ___ 

In contrast to the former resultal3'E2lO, the spin-gap 
phase spreads extensively toward the high-density region. 
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The spin-gap and phase-separation boundaries flow to- 
gether into the point J/t ~ 3.5 as n — > 1. We are not 
able to answer whether the spin gap survives in the n — > 1 
limit or not, because the numerical results become unsta- 
ble in the high density region where the phase boundary 
is close to the phase-separated state. 

In the low-density region, the phase boundary can be 
determined analytically by solving a two-electron prob- 
lem (see Appendix C). Then the asymptotic behavior of 
the phase boundary is obtained as 



2t/J c = cos(7m/2). 



(40) 



Note that the J c given by eq. ([l0|) in L — > oo limit is equiv- 
alent to the critical point where_the singlet pair forms a 
bound state in the ground stateEl This explains the fact 
that the spin-gap phase boundary overlaps the K p = 1 
line where the TL liquid behaves as free electrons, in the 
low-density limit. 

In order to check the consistency of our argument, we 
calculate the scaling dimensions for the singlet and the 
triplet excitations from eqs.(jl^) and (^). Then the av- 
erage of the renormalized scaling dimension ( p4| ) is taken 
so as to eliminate the logarithmic corrections as 



^.singlet _|_ ^triplet 



(41) 



x a and its finite-size effect are shown in Fig.|] and Fig|], 
respectively. The extrapolated data become 1/2 with 
error less than 0.2 %. 

In the spin-gap region (J > J c ), the asymptotic be- 
havior of the spin gap is obtained using the relation 
yo cx J — J c and eq. 



B. The t-J-J' Model 

Next, we analyze the ID t-J-J' model (||). According 
to Ref.|24|, the critical point for the J, J' — > limit is 
obtained by mapping the spin part onto the case of n = 1 , 
using the factorized wave function, 

J off = J(n l n. i+ i) S F + J'{n. t (l - ra i+ i)ra i+2 )sF, (42a) 
^cff = J'(n. 4 n i+ in i+2 )sF, ( 42b ) 

where (• • -)sf indicates the expectation value of the non- 
interacting spinless fermion. The effective ratio of the 
frustration a e s is then obtained as 



a e e(n, a) 



(1 + l/a)n 2 - s| - sf/a 
n 3 - (2sf + s|)n + 2s 1 2 s 2 



- 1 



(43) 



where s; = sin(^7rn)/Z7r. We can obtain the critical 
density n c where the spin gap vanishes, by compar- 
ing eq.(^3|) pjth the result of the frustrated spin chain: 
a c = 0.2411E1 For a = 1/2, we get n c = 0.7433. 

On the other hand, in the low-density limit, the criti- 
cal value for the spin-gap phase J c /t, can be analytically 



obtained by solving the two-electron problem as (see Ap- 
pendix C) 



4i/J c = 1 + 2a + VI + 4a 2 . 



(44) 



The meaning of this point is same as that of the t-J 
model. 

We show the phase diagram of a = a c case in Fig.||(b). 
The spin-gap phase boundary overlaps with the contour 
line of K p — 1 at almost all densities. This situation 
is quite resemble to that of the super-symmetpic t-J 
model with long-range hopping and interactionsE-3. In 
this model, there are no logarithmic corrections and the 
exact ground state is given by the Gutzwiller wave func- 
tion. This means that the charge degrees of freedom is 
free electrons (K p = 1), and the singlet and the triplet 
excitation spectra are degenerate for all densities. 

Figure ||(c) is the phase diagram at a = 1/2. The 
phase boundary starts from the critical value of the low- 
density limit (f44|), and bends at n ~ 2/3. It then flows 
into the critical point ( J/t, n) = (0, n c ). Thus the spin- 
gap phase with different origin (the Majumder-Ghosh- 
like dimer phase in the low-doping region, and the spin- 
gap phase in the large J/t region) have a single domain 
in the phase diagram. In contrast to the case of the t-J 
model (a = 0), the spin-gap phase boundary lies in the 
K p < 1 region so that there is no TS region in this case. 

Spin-gap phase may also exists for a > 1/2 cases. For 
example, in a — oo case, J' c /t — 1 in the dilute limit, and 
n c = 0.5752 in the J' —> limit. 

In spite of the deformation of the phase diagram, the 
critical value at the quarter-filling (n = 1/2) is almost 
independent of the strength of the frustration a, and is 
kept at J c /t ~ 2.7. Let us consider the reasoaior this us- 
ing an argument based on the g-ology modelQIS. In order 
to apply the g-ology model, we add the on-site Coulomb 
term to eq.(^|) and relax the constraint. The original 
Hamiltonian is restored when we set U — oo. Since the g- 
ology model is appropriate for the weak coupling case, we 
consider J' terms as corrections to the t-J model which 
belongs to the universality class of the TL model. Then 
their contributions to the ^-parameters, which are related 
to the spin-gap generation, are identified as 



Sgi± = Sg a = — J'(l + cos4/cf)- 



(45) 



For the quarter-filling, eq. (|45|) vanishes, so that the J' 
terms do not affect the renormalization flow of the spin 
part. Thus the frustration does not change the critical 
point at the quarter-filling within the scheme of the g- 
ology model. 



C. The Extended Hubbard Model 

The instability of the extended Hubbard model can 
be argued basedpCaa the g-ology model for weak coupling 
cases (U, V — * O)0'H. The g parameters which are used to 
determine the phase diagrams are identified as follows: 
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9i± = 9a = U + 2V A cos2fcp, 

g p = [7 + 2^(2 -cos 2fc F ) 



(46a) 
(46b) 



We have defined g v = gn — g 2 \\ T 92± where the up- 
per (lower) sign corresponds to v = p (y = a). Then, 
from the discussion given in Sec.|n|, the spin-gap phase 
boundary is determined by g±± = g a = 0. On the other 
hand, the contour line for K p = 1 is determined by the 
condition g p = 0, due to the following relations 



I 2ttv f +ff 4 || ±.g 4 _L +9u 
2irv F + 54|| ± g 4 ± - g v ' 



(47) 



where g 4 ± = U + 2V, g 4 \\ = 2V. 

Figure [?] shows the phase diagrams of the ID extended 
Hubbard model (||) for various electron fillings. They are 
obtained by analyzing the data of L = 12 systems. In 
the all cases, the slopes of the spin-gap phase boundaries 
and the K p = 1 contour lines near the origin of the U- 
V plain, are consistent with those predicted by the g- 
ology. For V < region, there is phase-separated state 
and its boundaries flow into (U,V) = (oo, — 2t) due to 
the equivalence of the XXZ spin chain in the large-?/ 
limitca. The spin-gap phase boundaries flow into these 
phase-separation boundaries. 

In n — 1/3 case, the spin-gap phase boundary and the 
K p = 1 contour line almost overlap near the solution of 
the two electron problem (see Appendix C), 



V r = - 



2U C 



U c /t- 



(48) 



This phenomenon is same as that of the low-density re- 
gion of the t-J(-J') model. 

At n = 1/2, the spin-gap phase boundary is clos e to 
U = 0, because the effect of V is canceled in eq.(46a). In 
this phase diagram, there are two regions with K p > 1. 
Besides of the spin-gap phase, a charge-gap phase exists 
for U, V > Q region due to the Umklapp scattering. Tip 
analysis for this instability will be reported elsewhereEl 
For.a.^ 2/3, a spin-gap phase appears in U/t,V/t > 



regionEfl. This is because the strong nearest-neighbor re- 
pulsion stabilizes the on-cite singlet pairs. The one of 
the striking feature in this phase diagram is that there 
are two phase separated states in the V/t ^> 1 region, 
and the spin-gap phase boundary flows between these 
two phase-separated states. In this region, the spin-gap 
phase boundary shifts to the large U side due to the 
strong finite-size effect. The phase-separated state in the 
U/t > side is considered as a mixture of Akp- and 2fcp- 
CDW phases. The stability of this phase is already ar- 
gued in Refjl^ by using the second-order perturbation 
theory. On the other hand, in the U/t < side, the sys- 
tem is separated into a 2fcp-CDW phase and a vacuum. 
These phase-separated states are illustrated in Fig.||. 

The consistency of the argument can also be checked 
as in the case of the t-J model. Fig shows the averaged 
scaling dimension (Ell) at n = 1/2 for V/t = 2 and 8 



cases, calculated by the data of L = 8, 12, 16 systems. 
Although the finite-size effect is large for the V ^> 1 
region, the extrapolated value become 1/2. 



V. CONCLUSION 

To conclude, we have studied critical properties of spin- 
gap phases in ID electron systems, considering the effect 
of the backward scattering in TL liquids by the renormal- 
ization group analysis. The phase boundary between TL 
liquids and spin-gap phases is shown to be determined by 
the singlet-triplet level crossing point. These excitation 
spectra are extracted by twisting boundary conditions, 
and identified by the discrete symmetries of wave func- 
tions. For this purpose, we have discussed symmetries of 
wave functions under parity and spin-reversal transfor- 
mations. We have applied the analysis to the extended 
Hubbard model and the t-J-(J') model, and clarified the 
spin-gap regions in the phase diagrams. The consistency 
of the our result has been checked by investigating the 
ratio of logarithmic corrections. Our results are also con- 
sistent with those of the g-ology model in the weak cou- 
pling limit, and of the two-electron problem in the dilute 
limit. 
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APPENDIX A: QUANTUM NUMBERS IN TWO 
NOTATIONS 

In the analysis of ID electron systems by Bethe-ansatz 
results with CFT, a different notation from Q'Tf ig often 
used to describe the quantum numbersOEflHo'EJ. In 
these notation the spin degrees of freedom is imposed 
only on down spins. In their definition, AN C is the change 
of the total number of electrons, and AN S is the change of 
the number of down spins. D c (D s ) denotes the number 
of particles moved from the left charge (spin) Fermi point 
to the right one. They are given by the eigen value of the 
number operator N r s as 



AA^ C = Y RT + Y LT 
AN s = N Rl +N Ll , 
2D C = Y RT - Y LT , 
2D C + 2D S = N Rl ~ N Ll . 



N- 



N) 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
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From eqs.(|12|) and flAl| ), the quantum numbers can be 
read as n p = AN c /2, n a = AN c /2 - AN S , m p = 2D C + 
D s , m a = — D a . One can also easily show the equivalence 
of the selection rule given by eq. ( |l3|) and the one written 
by this notationE3: 



AN C + AA S 
D c = (mod 1), 



D s = 



AN C , 



(mod 1). 



(A2a) 
(A2b) 



This relation is derived from the U — > oo limit of the 
Hubbard model. 



APPENDIX B: DERIVATION OF LOGARITHMIC 
CORRECTIONS 

Here we derive the logarithmic corrections given in 
eq.(|24|). Hereafter, we omit the spin index a. We con- 
sider perturbation terms which break the scale invariance 
as 



dr 
2^ 



XiOi(r). 



(Bl) 



Then the correction to the finite-size scaling is calculated 
within the first-order perturbation ascj 



Ei — Eq 



2nv 



2irv 



x i ~ ^ ' AjCji 



(B2) 



where \<pi) is the eigen state of Ei, and Cijk is a universal 
constant (OPE (operator product expansion) coefficient) 
fixed by a three-point function: 
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r 23 J r 31 



(B3) 



This coefficient can be derived from the following two 
ways. 

1. Abelian Bosonization 

The Lagrangian density of the spin part of eq. (j|) (the 
sine-Gordon model) is written as 



with 



1 2-kcP- 2na 2 



(B4) 

(B5a) 
(B5b) 



where the perturbation term Ci consists of the following 
two parts: one is a part of Gaussian model which denotes 
the deviation from the free case (K = 1). The other is the 
cosine term which stems from the backward scattering. 
They denote the effect of interaction between the left and 
the right Fermi points, and are written in the Euclidean 
space as 



Oo = -c?K- 1 [{v- x d T <j>f + (d^f 
O x = V^cosv^. 

Their coupling constants are given by 
2X = y (l), y/2X 1 =y 1 (l). 



(B6a) 
(B6b) 

(B7) 



For the SU(2) symmetric case yo{l) — yi(l), and yo(l) > 
0, the marginally irrelevant coupling is calculated from 
eq.(E3) as 



y In L + 1 



(B8) 



where the the bare coupling is defined as yo = J/o(0), and 
we have set I = InL. 

Now we consider the operators for the singlet and the 
triplet states as 



2 = V2 cos V2</>, 

3 = V2 sin V2(f>, 
e> 4 = exp(+iV2(9), 



(B9a) 
(B9b) 
(B9c) 



then the coefficients of their OPE with the marginal op- 
erators (|B6|) are obtained as 



C220 — C330 — - 
C221 = — C331 = 



K 
I 



a 



440 



f 

2K' 



C441 = 0. 



(BIO) 



Thus the scaling dimensions of the operators for si ngle t 
and triple t excitations are obtained from eqs.(B2), (jBq), 
and (B1C ). These are consistent with the results obtained 
by Gimarchi and SchulalJ. 



2. Non- Abelian Bosonization 

In the standard bosonization theory, systems are de- 
scribed in U(l) symmetric form, so that the explicit 
SU(2) symmetry in spin degrees of freedom is lost. In 
order to describe systems with higher symmetry, it is de- 
sirable to perform the calculation defining current fields 
that conserve the SU(2) symmetry. 

In SU(2) symmetric case, the system is described by 
chiral SU(2) currents that are defined as 



t ,f "a/3 1 



(Bll) 







where cr = [a 1 , a 2 , a 3 ] are the Pauli matrices and r — 
(R, L). The chiral SU(2) current Jr has a conformal 
dimension (A+,A~) = (1,0) and J L has (A+,A _ ) = 
(0,1). The three components of J r obey commutation 
relations known as the Kac-Moody algebra with central 
charge k: 



fc/2 



i£ijldJ l r {w) 



(z - w) 2 



2 « 



reg., (B12) 



where em is the anti-symmetric structure factor. For 
spin-s systems, there is a relation k = s/2. If a sys- 
tem is described by this current algebra, the system 
belongs to the universality class of the Wess-Zumino- 
NovikoAt-iW|itten non-linear a model with topological cou- 
pling /£3'E3. 

In this case, the scaling dimension is 



2s r (s r + 1) 
2 + k ' 



(B13) 



where s r = 0, 1/2, fc/2. Therefore, the lowest energy 
spectra for the singlet and the triplet excitations are 



singlet 



X 



triplet 



(B14) 



Now let us consider the correction in the presence of a 
marginal operator (x = 2)EH which is given by 



= J L -J v 



(B15) 



The marginal operator O is proportional to Sx ■ S-r where 
S r is the SU(2) charge, and S = Sl + Sr, is the spin of 
the state fa. Letting the degrees of S and S r are s and 
s r , respectively, the expectation value becomes 

(fa\S L ■ S R \fa) = ~(s(s + 1) - s L (s L + 1) - s R (s R + 1)). 

(B16) 

Here, sl = sr = 1/2 and s — for the singlet and 
s = 1 for the triplet. Thus the ratio of the logarithmic 
corrections is calculated as 3 : —1. 



APPENDIX C: DILUTE LIMIT 

In the low-density limit, a many-body problem may 
be reduced to a two-body problem. Here we consider a 
critical point where a bound electron pair become stable 
in the ground state, and a singlet-triplet level-crossing 
takes place. We perform the calculation following thp 
approach of H. Q. Lin, which was used for the 2D caseEi 

In order to take the constraint of the t-J(-J') model 
into account, we relax the restriction, and add the on-site 
Coulomb term as 



n = n 



(CI) 



The result of the original Hamiltonian can be obtained 
when we set U = oo in the end of the calculation. 

It is well known for a two-body problem that the 
ground state is a singlet p^fer as the bottom of the energy 
band has no degeneracyE3C3. This is consistent with the 
argument in Sec. III. The wave function in this system 
can be written using the basis A as 



(C2) 



where f(i,j) = f(j,i) for the singlet (T = 1) as shown 
in Sec. III. The Schrodinger equation for the singlet wave 
function is 

Ef(i,j) = $>;/a,j)+ W(i,0] + PSij - Jij]f(i,j), 
i 

(C3) 

with tij — — t<5ij_jM and = J8u-j\ x+otJ5\i-j\2 where 
a denotes the strength of the frustration a = J' /J. The 
Fourier transformation of eq. flC3j ) is given by 



k), (C4) 



Ef(k 1 ,k 2 ) - [t(k 1 ) + t(k 2 )}f(k 1 ,k 2 ) 

Lj2[U-J(k)]f(h + k,k 2 



where 



/(fci,fc 2 ) = ^/(i,i)e- ife ^- ifc ^ i 

i'j 

t(k) = — 2icosfc, 
J(k) — 2 J(cos k + a cos 2k) . 



(C5) 

(C6) 
(C7) 



Next, we introduce center of mass and relative mo- 
menta by Q = ki + k 2 , q = (k\ — k 2 )/2, and redefine the 
functions as 



f Q (k) = f(k u k 2 ), 
e Q (q) = t(Q/2 + q)+t{Q/2-q). 

Then we get 

, , , fEJgW4E t J(g-W) 
jQii 5 r\ • 

where 

J(q — k) = 2 J(cos q cos k + a cos 2q cos 2k). (Cll) 
Note that the terms that contain sin are omitted in 



(C8) 
(C9) 



(CIO) 



eq.( |CllD , because they give no contribution due to their 
symmetry. Now we define the following variables, and 
iterate them as 

= UIq.qCq - 2JUI lfi d - 2aJUI oa C 2 , (C12a) 
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d = ^-^2f Q (q)cosq 

= h.oCo - 2JI 2 ,qC 1 - 2aJIi.iC 2 
Ci = — ^2 f Q (q) cos 2q 



where 



IoaCo ~ 2JI\ i \C\ — 2aJI§,iCi, 

1 ^ cosm 9 cos " 2 9 
= IZ^ E-e Q (q) ■ 



(C12b) 



(C12c) 



(C13) 



The criterion that eq. (|C12|) have a solution is 
det 



1 - Uhfi 2JUh fi 2aJUI ,i 
-h,o 2JI 2 , + 1 2aJI h i 
-J ,i 2JI X i 2aJI Q:2 + 1 



0, (C14) 



where all I m , n can be related to /o,o as follows, 

h.o = (1 - EI 0fi )/U, 

h,o = -Ehjo/U, 

7 3)0 - (1 - 2EI 2<0 )/8t, 

h,o = -Eh, /U, (C15) 

/l,l = 2/3^0 — /l,0i 
/o,l = 2/2,0 — Iq,0i 

I0.2 = Iq,o — 4/2,0 + 4/4,0, 

and /o,o diverges. Then setting U — 00, we get the rela- 
tion between the singlet-state energy and the parameters 
of the model as 

At/ J = ~z(Aaz 2 -2a+l) 



+ y/z 2 {Aaz 2 -2a + l) 2 - 4a(2z 2 - 1), (C16) 

where z = E/At. For the singlet pair with Q = 0, the 
energy is given by E = — At + B where B is the binding 
energy. At the critical point where the singlet pair be- 
comes stable, the binding energy becomes B = 0. Then 
we get the solution ( p4| ) without size dependence. In 
a = case, we obtain J c = 2t. 

In the case of the extended Hubbard model, the so- 
luti on c an be obtained by setting (J, a) = (— V, 0) in 
eq.( C14 ), and leaving U finite. The result is 



V = 



2U 



z{U/t-Az)' 



(C17) 



For U — > limit, it coincides with the spin-gap phase 
boundary and the K p = 1 contour line predicted by the 
g-ology: V = -U/2. 

Finally, we consider the singlet-triplet level-crossing 
point in the dilute limit. In the system with anti-periodic 
boundary conditions, the bottom of the energy band is 
degenerate, so that a level crossing may take place. For 
the triplet state, the last term of eq.(C4) vanishes due 



to the symmetry of the wave function: f(i,j) — —f(j,i) 
(T = —1). Therefore, the triplet state is always non- 
interacting. This means th at th e level-crossing point can 
be obtained as a solution (C16) for E — £q = q(it / L). In 
this case, the density dependence of the critical point of 
the t-J-J' model can be expanded as 



J c (n) = J c (0) + A{a)n 2 + C(n 4 ), 



(C18) 



where J c (0) is same as the solution for the ground 
state. Therefore, the spin-gap phase boundary in the 
low-density limit coincides with the critical point for the 
bound electron pair in the ground state, and its curve is 
the square-root type in the J/t-n plane. For a = case, 
we obtain eq.(|40|). These solutions reflect the shape of 
the band structure. 
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FIG. 1. Spectral flow of the ID t-J model at J/t = 2 with 
length L = 8 and electron number N — 4. These energy spec- 
tra are the lowest two levels for each wave number k — 2nI/L. 
The marked spectra correspond to the physical states written 
in this figure. 
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FIG. 2. Singlet and triplet excitation energies of the ID t-J 
model for L — 16 system at n — 1 /2. These excitation spectra 
can be identified by symmetries of their wave functions. 



12 




0.015 



FIG. 3. Size dependence of J c /t determined by the intersec- 
tions of the excitation spectra for L = 8, 12, 16, 20 systems at 
n= 1/2. These points are fitted by the form A+B/L 2 +C/L 4 . 
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FIG. 4. Extrapolated value of (x s a inglct +3x t J iplct ) /4 and the 
scaling dimensions for the singlet and the triplet excitations 
for L — 16 system at n — 1/2. 
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FIG. 6. Phase diagrams of the ID t-J-J' model at (a) 
a = 0, (b) a = q c , (c) a = 1/2 (TL: TL phase, SG: spin-gap 
phase, PS: phase-separated state). In the spin-gap phase 
where the backward scattering is attractive, the singlet ex- 
citation becomes lower than the triplet (see FIG-CM). The 
contour lines of K p are calculated by the data of L = 16 
system. 
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FIG. 7. Phase diagrams of the ID extended Hubbard 
model determined by the data of L = 12 systems at (a) 
n = 1/3, (b) n = 1/2, (c) ra = 2/3 (SDW (TS): TL liquid 
phase with K p < 1 (A" p > 1), CDW (SS): spin-gap phase 
with 7f p < 1 (K p > 1), PS: phase-separated state). 
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FIG. 8. Two phase-separated states which appear in the 
V/t > 1 region of Fig.|(c). (a) is located in U/t > side, (b) 
lies in U/t < side. Spin-gap phase boundary exists between 
these two phase-separated states in the phase diagram. 
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FIG. 9. Extrapolated value of (s^ inglct + 3a:^ iplct )/4 and the 
scaling dimensions for the singlet and the triplet excitations 
for L — 16 system at n — 1/2. (a) is the case of V/t = 2, and 
(b) is the case of V/t = 8. 
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TABLE I. Discrete symmetries of wave functions for dif- 
ferent two bases (V: space inversion, T: spin reversal, fc:wave 
number BC: boundary conditions). The upper (lower) sign 
denotes the case of N/2 =even (odd). These correspon- 
dences are explained by the Perron-Frobenius theorem and 
the bosonization theory. 
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